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Abstract 

Knowledge of the B s decay fraction of the T(5S) resonance, fs, is important for B s meson studies at 
the T(55) energy. Using a data sample collected by the CLEO III detector at CESR consisting of 0.423 
fb _1 on the T(5S) resonance, 6.34 fb _1 on the T(4S) and 2.32 fb _1 in the continuum below the T(4S), 
we measure B(T(5S) <f>X) = (13.8±0.7l?;|)% and B(T(4S) -»• <f>X) = (7.1±0.1±0.6)%; the ratio of 
the two rates is (1.9±0.1^o !)■ This is the first measurement of the <f> meson yield from the T(5S). Using 
these rates, and a model dependent estimate of B(B S — > (f>X), we determine fs = (24.6 ± 2.9~^_ 1 l'^)%. 
We also update our previous independent measurement of fs made using the inclusive D s yields to 
now be (16.8±2.6l" 3 i ' 4 ')%, due to a better estimate of the number of hadronic events. We also report the 
total T(55) hadronic cross section above continuum to be a(e + e~ — > T(55)) = (0.301 ±0.002 ±0.039) 
nb. This allows us to extract the fraction of B mesons as (58. 9± 10.0±9.2)%, equal to 1-fs- Averaging 
the three methods gives a model dependent result of fs = (2ll|)%. 

PACS numbers: 13.20.Hc 
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I. INTRODUCTION 



The putative T(5S) resonance was discovered at CESR long ago by the CLEO [1] and 
CUSB [2] collaborations by observing an enhancement in the total e + e~ annihilation cross- 
section into hadrons at a center-of-mass energy of about 40 MeV above the B*B S production 
threshold. Its mass and its production cross-section were measured as (10.865±0.008) GeV/c 2 
and about 0.35 nb, respectively [1]. 

The possible final states of the T(5S) resonance decays are: B*B S *, B S B S *, B S B S , B*B*. 
BB , BB, BBir, BB tt, B*B tt, BBitit, where the B and 7r mesons can be either neutral or 
charged. Several models involving coupled channel calculations have predicted the cross-section 
and final state composition [3, 4]. The Unitarized Quark Model [4], for example, predicts that 
the total bb cross-section at the T(5S) energy is dominated by B*B and B*B S production, 
with B s production accounting for about 1/3 of the total rate. The original ~116 pb" 1 of data 
collected at CESR did not reveal if any B s mesons were produced. By measuring the inclusive 
D s production rate using CLEO III data we previously measured the fraction of B s meson 
production at the T(5S), fs, to be (16.0 ±2.6 ±5.8)% of the total bb rate [5]; this measurement 
uses a theoretical estimate of B(B S — > D S X) = (92 ± 11)% [5]. 

B s production was confirmed in a second CLEO analysis [6] that fully reconstructed B s 
meson decays. In addition, this analysis showed that the final states are dominated by the 
B*B S decay channel. These results have been confirmed by the Belle collaboration [7]. A third 
CLEO analysis of the exclusive B reconstruction at the T(5S) showed that the BBX final 
states in T(5S) decays are dominated by B*B with a considerable contribution from BB and 
B*B final states [8]. 

Knowledge of the B s production rate at the T(5S) resonance, fs, is necessary to compare 
with predictions of theoretical models of 6-hadron production. More importantly, fs is also 
essential for evaluating the possibility of B s studies at the T(5S) using current 5-factories, and 
for future e + e~ Super- B Factories, should they come to fruition, where precision measurements 
of B s decays are an important goal [9] . However, the determination of fs in a model independent 
manner requires several tens of fb _1 [10]. In this paper, we improve our knowledge of the B s 
production at the T(5S) by using inclusive meson yields to make a second model dependent 
measurement of B(T(5S) — > B*B S *). 

We choose to examine meson yields because they will be produced much more often in B s 
decays than in B decays. The branching fractions of B s to D s mesons and B to D mesons are 
of order 1, while the rates of B s to D mesons and B to D s mesons are on the order of 1/10. 
We also know that the production rate of mesons is one order of magnitude higher in D s 
decays than in D° and D + decays as measured using recent CLEO-c data [11]. These results 
are shown in Table I. The rate of D mesons into mesons is only 1%, while the rate of D s 
mesons into mesons is 16%. 

Since the B — > <pX branching ratio already has been measured to be (3.42 ± 0.13)% [12], 
we expect a large difference between the yields at the T(5S) and at the T(4S), due to the 

presence of B s , that we will use to measure the size of the B^B S ^ component at the T(5S). 
The analysis technique used here is similar to the one used in [5]. We also present an update to 
the first measurement of fs that used D s yields [5], and we report on the measurement of the 
total T(5S) hadronic cross-section. When we discuss the T(5S) here, we mean any production 
above what is expected from continuum production of quarks lighter than the b at an e + e~ 
center-of-mass energy of 10.865 GeV. 
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TABLE I: D°, D + and inclusive branching ratios into cp mesons. The D meson decays used 281 
pb _1 of data at the 0(3770), while the D s decays were measured using 195 pb _1 at or near 4170 MeV. 

D° 1.05 ±0.08 ±0.07 
D + 1.03 ±0.10 ±0.07 
D+ 16.1 ±1.2 ±1.1 



The CLEO III detector is equipped to measure the momenta and directions of charged 
particles, identify charged hadrons, detect photons, and determine with good precision their 
directions and energies. It has been described in detail previously [13, 14]. 

II. MEASUREMENT OF B(T(5S) B^B^^) USING MESON YIELDS 

A. Data Sample and Signal Selection 

We use 6.34 fb -1 integrated luminosity of data collected on the T(AS) resonance peak and 
0.423 fb -1 of data collected on the T(5S') resonance (E C m = 10.868 GeV). A third data sample 
of 2.32 fb -1 collected in the continuum 40 MeV in center-of-mass energy below the T(4S') is 
used to subtract the four-flavor (u, d, s and c quark) continuum events. 

Hadronic events are selected using criteria based on the number of charged tracks and the 
amount of energy deposited in the electromagnetic calorimeter. To select "spherical" 6-quark 
events we require that the Fox- Wolfram shape parameter [15], R2, be less than 0.25. meson 
candidates are looked for through the reconstruction of a pair of oppositely charged tracks 
identified as kaons. These tracks are required to originate from the main interaction point and 
have a minimum of half of the maximal number of hits in the tracking chambers. They also must 
satisfy kaon identification criteria that uses information from both the Ring Imaging Cherenkov 
(RICH) and the ionization loss in the drift chamber, dE/dx, of the CLEO III detector. Kaon 
identification has been described in detail previously [5]. 

B. Meson Yields From T(5S) and T(4S) Decays 

All pairs of oppositely charged kaon candidates were examined for candidates if their 
summed momenta is less than half of the beam energy. Instead of momentum we choose to 
work with the variable x which is the momentum divided by the beam energy, to remove 
differences caused by the the change of energies between continuum data taken just below the 
T(4S), at the T(4S) and at the T(5S). The K + K~ invariant mass distributions for x < 0.5 are 
shown in Fig. 1. 

1. K + K~ Invariant Mass Spectra and Yields 

A crucial aspect in this analysis is to accurately model the signal and background shapes. 
For the signal we use a Breit-Wigner signal shape convoluted with a Gaussian. The width of 
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FIG. 1: The invariant mass distributions of the candidates with x < 0.5 from: (a) the T(4S) on- 
resonance data (b) the continuum below the T(4S) resonance data, and (c) the T(5S) on-resonance 
data. The solid line is the fit to the signal and background shapes explained in the text. 



the Breit-Wigner function was fixed to the natural width of the <fi meson, I\ = 4.26 MeV [12]. 
It is convoluted here with a Gaussian function to allow the integration of the detector resolution 
into the signal function. A second Gaussian is added for an adequate fitting of the tails. The 
form of the signal fitting function, where the dependent variable ( is the invariant mass of our 
K + K~ candidates, is given by 

s(0 = HO + [ C+5a f(y)g((-y)dy. (l) 
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Here, at every physical point £, the Breit-Wigner function 

p / \ i a 1 -t 1 




(2) 
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The integration over y is between the limits C — 5a, and C + 5cr. Ax and y are the area and 
mean value of the Breit-Wigner, and 02 is the standard deviation of the Gaussian. h(Q is a 
second Gaussian added to fit the tails. 

The background shape is given by a function, chosen to model the threshold, that has the 
following functional form: 



where A is the normalization, Co is the threshold, p is the power of a polynomial about the 
turn-on point, and c\ and C2 are linear and quadratic coefficients in the exponential. Because 
this function has a sharp rise followed by an exponential tail, we are able to accurately describe 
the threshold behavior at low invariant mass close to the kinematic limit. 

We show the invariant mass of the K + K~ candidates in 9 different x intervals (from 0.05 to 
0.50) for all the data samples in Fig. 2, Fig. 3 and Fig. 4. We determine the raw yield of 0's 
shown in Fig. 5 (a), (b) and (c) from fits of the invariant mass distributions to the functions 
S(Q and -B(C) defined above. All the parameters in S(() in all the fits, except the yields, were 
fixed to the values obtained when fitting the corresponding distributions from the sum of the 
three data samples collected at the three different beam energies in each x interval separately. 
The parameters describing B((), on the other hand, are allowed to float except for Co which is 
fixed to twice the K + mass. The yields are listed in the second and third columns of Table II 
and Table III. 

The systematic errors on these yields are estimated to be ±4%, and are obtained by varying 
the fitting techniques. One variation is to float the fitting parameters at the different beam 
energies, instead of fixing them. We also used different functions including a Breit-Wigner 
signal shape with and without convolution, or a double Guassian shape for the signal, and for 
background either other threshold functions or polynomials. 

2. Continuum Subtraction 

The numbers of hadronic events and <fi candidates from the T(5S) and the T(4S) resonance 
decays are determined by subtracting the scaled four-flavor {u, d, s and c quarks) continuum 
events from the T(4S) and the T(5S) data. The scale factors, S n s, are determined using the 
same technique described in [5], where 



with L n s, L cont , E n s and E cont being the collected luminosities and the center-of-mass energies 
at the T(nS) and at the continuum below the T(4S). We find: 



B(Q = A-((- (of ■ exp [ Cl (C - Co) + c 2 (C - Co) 2 ], 



(4) 
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FIG. 2: The K + K~ mass combinations from T(4S) on-resonance data, fitted to the sum of a Breit- 
Wigner signal shape centered at the nominal 4> mass, convoluted with a Gaussian to describe the 
detector resolution, and a second Gaussian distribution for a better parametrization of the tails. The 
background is parameterized by a threshold function (see text). These distributions are in the x 
intervals: (a) 0.05 < x < 0.10, (b) 0.10 < x < 0.15, (c) 0.15 < x < 0.20, (d) 0.20 < x < 0.25, (e) 
0.25 < x < 0.30, (f) 0.30 < x < 0.35, (g) 0.35 < x < 0.40, (h) 0.40 < x < 0.45, (i) 0.45 < x < 0.50. 
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FIG. 3: The K + K~ mass combinations from the continuum below the T(4S), fitted to the sum of a 
Breit-Wigner signal shape centered at the nominal (j) mass, convoluted with a Gaussian to describe 
the detector resolution, and a second Gaussian distribution for a better parametrization of the tails. 
The background is parameterized by a threshold function (see text). These distributions are in the 
x intervals: (a) 0.05 < x < 0.10, (b) 0.10 < x < 0.15, (c) 0.15 < x < 0.20, (d) 0.20 < x < 0.25, (e) 
0.25 < x < 0.30, (f) 0.30 < x < 0.35, (g) 0.35 < x < 0.40, (h) 0.40 < x < 0.45, (i) 0.45 < x < 0.50. 
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FIG. 4: The K + K~ mass combinations from the T(5S), fitted to the sum of a Breit-Wigner signal 
shape centered at the nominal (f> mass, convoluted with a Gaussian to describe the detector resolution, 
and a second Gaussian distribution for a better parametrization of the tails. The background is 
parameterized by a threshold function (see text). These distributions are in the x intervals: (a) 
0.05 < x < 0.10, (b) 0.10 < x < 0.15, (c) 0.15 < x < 0.20, (d) 0.20 < x < 0.25, (e) 0.25 < x < 0.30, 
(f) 0.30 < x < 0.35, (g) 0.35 < x < 0.40, (h) 0.40 < x < 0.45, (i) 0.45 < x < 0.50. 
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FIG. 5: Raw <p yields from: (a) the T(4S) data, (b) the continuum below the T(4S) data, (c) the 
T(5S) data. 



S 4S = 2.713 ± 0.001 ± 0.027 



(6) 



and 



S 55 = (17.18 ±0.01 ±0.17) • 10~ 2 . (7) 
We estimate the systematic error on these scale factors by obtaining them in a different 
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TABLE II: (ft yields from the T(4S) (N^ 4S)on ), from the continuum below the T(4S) (N l cont ), and 
from the continuum subtracted T(4S) (^Vx(4S))- ^ so listed are the (ft reconstruction efficiencies (e l 



and the partial Y(4S) 
yields are ±4%. 



4>X branching fractions as a function of x. The systematic errors on the 



■) 



A^ 



Y(45)on 



AT 



cont 



iV T(45) 



0.05-0.10 
0.10-0.15 
0.15-0.20 
0.20-0.25 
0.25-0.30 
0.30-0.35 
0.35-0.40 
0.40-0.45 
0.45-0.50 



4938.1 ± 135.2 
22848.9 ± 228.7 
33207.4 ±269.1 
26020.4 ± 237.6 
15364.1 ± 178.8 

8107.8 ± 128.2 
4263.5 ± 95.9 
2475.3 ±75.1 

1256.9 ±57.6 



462.4 ±41.0 
2179.0 ±69.0 
3608.7 ±84.1 
3218.7 ±80.7 
2374.2 ± 69.4 
1597.6 ± 57.8 
1113.4 ±48.6 
786.6 ±42.3 
443.6 ± 30.9 



3683.5 ± 175.1 8.0 

16937.3 ± 295.5 31.9 

23417.2 ± 352.8 47.1 

17288.3 ±323.1 51.3 
8923.0 ±259.7 49.1 
3773.5 ± 202.5 42.3 
1242.9 ± 162.9 35.8 
341.2 ±137.1 21.6 
53.5 ±101.8 13.4 



1.46 ± 0.09 
1.68 ± 0.04 
1.57 ±0.03 
1.06 ± 0.03 
0.58 ± 0.02 
0.28 ± 0.02 
0.11 ±0.01 
0.05 ± 0.02 
0.01 ± 0.02 



TABLE III: (ft yields from the T(5S) (A^ (5S)o J, from the continuum below the T(4S) (N l cont are 
the same as in the previous table), and from the continuum subtracted T(5S) continuum subtracted 
(A^ 5 ^). Also listed are the (ft reconstruction efficiencies (e* which are taken to be the same as in the 
previous table), and the partial T(5S I ) — ► (ftX branching fractions as a function of x. The systematic 
errors on the yields are ±4%. 

^£r) N ^S)on ^cont ^ (5S ) g%) 

0.05-0.10 214.7 ±26.2 462.4 ±41.0 135.3 ±27.1 8.0 2.7 ± 0.6 

0.10-0.15 1034.8 ±46.0 2179.0 ± 69.0 660.4 ±47.5 31.9 3.3 ± 0.2 

0.15-0.20 1432.5 ± 53.0 3608.7 ±84.1 812.4 ± 54.9 47.1 2.8 ± 0.2 

0.20-0.25 1151.9 ±47.3 3218.7 ± 80.7 598.9 ±49.3 51.3 1.9 ± 0.2 

0.25-0.30 789.2 ±38.6 2374.2 ± 69.4 390.2 ±40.4 49.1 1.3 ±0.1 

0.30-0.35 382.3 ±27.9 1597.6 ± 57.8 107.8 ± 29.7 42.3 0.4 ±0.1 

0.35-0.40 261.2 ±22.7 1113.4 ±48.6 70.0 ±24.1 35.8 0.3 ±0.1 

0.40-0.45 169.7 ± 19.1 786.6 ± 42.3 34.5 ± 20.4 21.6 0.3 ± 0.2 

0.45-0.50 75.0 ± 13.6 443.6 ± 30.9 -1.2 ±14.6 13.4 -0.1 ±0.2 



manner. Here we measure the ratio of the number of charged tracks at the different beam 
energies in the 0.6 < x < 0.8 interval. The lower limit on the x interval used here is determined 
by the maximum value tracks from BB events can have, including smearing because of the 
measuring resolution, and the upper limit is chosen to eliminate radiative electromagnetic 
processes. Since the tracks should be produced only from continuum events, we suppress 
beam-gas and beam-wall interactions, photon pair and r pair events using strict cuts on track 
multiplicities, event energies and event shapes. Since particle production may be larger at the 
higher T(5S f ) energy than the continuum below the T(AS), we apply a small multiplicative 
correction of (1.016±0.011)%, as determined by Monte Carlo simulation to the relative track 
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yields. Track counting gives a 1% lower value for S$s and we use this difference as our estimate 
of the systematic error. Note, that the error on the beam energy (0.1%) has a negligible effect. 

These numbers are updated from those reported in [5], due to an increase in Monte Carlo 
statistics, and the use of a more precise calibration of the beam energy. 



3. 4> Reconstruction Efficiency 



The x dependent (ft detection efficiency shown in Fig. 6 is determined by reconstructing and 
fitting (ft candidates in more than 8 million simulated BB generic decays. The low reconstruction 
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FIG. 6: (ft reconstruction efficiency from more than 8 million BB Monte Carlo simulated events. The 
errors are statistical only. 



efficiency at small x is due to the fact that as the (ft becomes less energetic, it becomes more 
probable that it is formed of a slow kaon (with momentum below 0.2 GeV/c) whose detection 
is very inefficient, since it is likely to be absorbed in the beam-pipe or vertex detector material. 
The efficiency for large x (0.4 < x < 0.5) is somewhat lower than naively expected because of 
the i?2 cut of 0.25. Due to the low efficiency and large backgrounds in the first bin, x < 0.05, 
we do not measure the (ft yield in this interval, but will rely on a model of (ft production to 
extract a value. 
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4- 4> Branching Ratios 



To find the number of hadronic decays produced above the four-flavor continuum, denoted 
as N^(nS)i we multiply the number of events found in the continuum below the T(4S), N^^, 
by the S^s an d Sss scale factors, and subtract them from the number of hadronic events found 
at each resonance, N^ nS y. 

jV* e | 5) = 7V™ 4S) - S 4S • Nf (4S) = (6.42 ± 0.01 ± 0.26) • 10 6 (8) 

N*H S) = N^ 5S) - S 5S ■ Nf (5S) = (0.127 ± 0.001 ± 0.016) • 10 6 . (9) 

Using these numbers together with N^/^ (iVr(5S)) anc ^ e% which are the yield and 
reconstruction efficiency in the i-th x interval respectively, we measure the inclusive partial 
decay rate of the T(nS) resonance into mesons in the i-th x interval as follows: 

B *' T <" S » - « = N^ S) x K+K-) fV) (10) 

The results are listed in the last column of Table II and Table III and shown in Fig. 7. 
Summing the results for 0.05 < x < 0.50 we measure 

B^ X>0 - 05 \T(AS) -> (f>X) = (6.82 ± 0.09 ± 0.55)%. (11) 

As the rate in the first interval 0.00< x <0.05 is experimentally not accessible in this analysis, 
we estimate the branching fraction in this interval using our Monte Carlo simulation to model 
4> production. This interval has 3.6% of the total yield, which corresponds to a branching 
fraction of 0.24%. We take the error as equal to the value. The total production rate then is 

B(T(4S) -> <j>X) = (7.06 ± 0.09 ± 0.60)%, (12) 
that when divided by two gives the B meson branching fraction into mesons of 

B(B -> <j>X) = (3.53 ±0.05 ±0.30)% , (13) 

which is in good agreement with the PDG [12] value of (3.42 ± 0.13)%. For the T(5S) we 
measure for x > 0.05 

B (x>om \T(5S) -> 0X) = (12.9 ± 0.7l?i)%. (14) 
In addition, we find for x > 0.05 

B(T(5S) <j>X) | (15) 



B(T(4S) -> <f>X) 

This result is more than 9.6<r significant, including both the statistical and the systematic 
errors, thus demonstrating an almost factor of two larger production of mesons at the T(5S) 
than at the T(4S). 
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FIG. 7: The efficiency corrected <f> branching ratio in each x interval from the T(5S) resonance decays 
(triangles), and from the T(4S) resonance decays (open circles). Errors are statistical only, except for 
the first bin where the T(4S) data point is a theoretical estimate based on Monte Carlo, with an error 
equal to its value. 

5. Determination of The Statistical and Systematic Uncertainties 



Since the measurements presented in the previous section depend on a large number of 
parameters, the corresponding errors on those measurements can have highly correlated con- 
tributions. In this section, we explain the method we used to extract the statistical and the 
systematic errors on our measurements by taking into account the sources of correlations be- 



14 



tween the different parameters. 

We measure the production rate at the T(4S) and T(5S) using 

1 10 /V* 

^>0.05) (T(nS) _ m = E 

iV T(nS) ' °1 i=2 C 

1 10 jV' — Q . M i 

I -sr^ JV T(nS)on D « JV cont 



l (na)on " ^'cont (16) 



( N T(nS)<m ~~ ^ ' ^cont) ' #1 i=2 

The quantities in this equation are defined as: 

• ^rfnS) * s the number of hadronic events at the T(nS) energy after continuum subtraction. 

• N^^ nS ^ on and N^ ont are respectively the number of hadronic events from the data taken 
at the T(nS) energy and from the continuum data taken at 30 MeV below the T(4S) 
resonance. 

• B 1 is the B((f) -> K + K-) branching fraction of 49.1% [12]. 

• -^r(nS) i s the number of candidates in the i-th x interval from the T(nS) resonance 
events (i.e. after continuum subtraction). 

• Ny./ nS \ on and A^ ont are the number of candidates in the i-th x interval from the T(nS)- 
ON resonance data and from the data taken at the continuum below the T(4S) respec- 
tively. 

• e l is the reconstruction efficiency of candidates in the i-th x interval. 

• S n are the continuum subtraction scale factors described in the previous sections. 

The ratio of the production rate from the T(5S) to the production rate from the T(4S) 
is given by 

1 V-IO N T(5S) 

_ g(»>o-Q5)(T(5S) -> 0X) — 

B(->o.o5)( T (4S) ^0X) i v io ^T ( 45) 1 ' 

1 "T^IO jV T(5g) n~ ,g5 * jV cont 

= (^( 5S)o n-g5*^ nt )-Bl ^ = 2 

1 y^lO jV Y(45) n~ S4 * jV cont 

K(4 S )o„- S 4*^o nt )^l ^ = 2 

The independent parameters are: S4 and S 5 , e\ Ny, nS , and A^ ont . The error on e l includes, as 
shown in Table IV, a 2% error on the tracking efficiency, a 2% error on the particle identification, 
both per track. In addition, it has a contribution of 1% from the error on the fitting method. 
The total systematic error on the reconstruction efficiency is therefore equal to 5.7%, and is 
correlated among all the bins. 

We determine the errors, including all the correlations, by using a Monte Carlo method 
where we generate each independent quantity as a Gaussian distribution using the estimated 
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TABLE IV: Systematic errors on the <f> detection efficiency. 



Systematic error 



Track finding (per track) 2 
Particle identification (per track) 2 
Fitting techniques 1 



Total 



5.7 



mean as the central value and the estimated error as the width. Then we evaluate the relevant 
measured quantities. Applying this method to R gives the spectra shown in Fig. 8. We do 
this separately for the statistical and systematic errors and then combine them for the total 
uncertainty. This is the probability distribution for this ratio. It is non-Gaussian. 



80000 — 



1630706-020 

r 




FIG. 8: Probability distributions for R smeared by its: (a) statistical error (solid line), (b) systematic 
error (filled triangles) and, (c) total error (filled squares). 



We extract the statistical, systematic and total errors on our measurements from the values 
of these distributions corresponding to 68.3% (one standard deviation) of the area under the 
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corresponding curves (either statistical, systematic or total) above and below the maximum 
position. We measure 

R = 1.9 iO.llgj . (18) 
C. Measurement of the Y(5S) — ► B^B^ Branching Fraction 

Eq. (15) and the spectrum in Fig. 7 demonstrate a significant excess of mesons at the T(5S). 

From these results, we can determine the fraction of the T(5S) that decays into B^fi'*', which 
we denote as fs, in a model dependent manner. We assume that the yields at the T(5S) 
comes from two sources, B and B s mesons. The equation linking them is 

B(T(5S) - 4>X)/2 = f s ■ B(B S - <t>X) + (1 • B{T(AS) - 0X). (19) 

We restrict our consideration to the interval 0.05 < x < 0.50. We obtain fs via the equation 
[16] 

_ jg(*>o-05)( T (5,g) -> 0X) - B^ >om \T(4S) -> 0X) 
/s ~ 2i3^>°- 05 )(5 s 0X) - B(*> - 05 \T{4S) 0X) ' ^ ^ 

The branching fractions i3( a;>a05 )(T(5 1 S) -> 0X) and i3( :c>a05 )(T(45) -> 0X) are given in Eqs. 
11 and 14. It is necessary to estimate B^ X>0 ' 05 \B S — + 0X). We first show that most of the 0's 
in b decay result from the decay chain B — > (D or D S )X, D or D s — > 0X. The decay rates for 
£> — > are tabulated by the PDG [12]. CLEO-c has measured, in a companion paper [11], 
the branching ratios for D°, D + and Df mesons into mesons. These results are listed in the 
second line of Table V. 



TABLE V: B decay inclusive branching ratios into D°, D + and Df mesons, the inclusive D decay 
branching ratios into (f> mesons and the resulting B — > <pX product branching ratios. 





D° (%) 




- (%) 


D+ (%) 


B(B -► (D or D S )X) 


64.0 ±3.0 


221 


!± 1.4 


8.6 ± 1.2 


B((D or D s ) 0X) 


1.05 ±0.10 


1.03 


±0.12 


16.1 ± 1.6 


B(B -► (D or D S )X) • or 0X) 


0.67 ± 0.08 


0.23 


±0.03 


1.4 ±0.2 



The rates on the last line of Table V give the resulting yields from the decay of B into 
each of the individual charm mesons and the subsequent decay of the charm mesons into a <fi. 
The sum of these product rates is (2.3±0.3)%. The difference with the measured B — > <pX 
rate given above in Eq. (13) is an additional unaccounted for (1.2±0.4)%. This rate can be 
accounted for from the production of charm baryons or merely by fragmentation processes. The 
sum of B — > DX and D S X branching ratios Table V is (95.4 ± 3.5)%. We now assume that 
the number of D plus mesons produced in B s decays is the same as in B decays. Using 
our previous estimate of B(B S — > D S X) = (92 ± 11)% [5], we have (8±11)% of the B s rate into 
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charmed mesons that must be accounted for by a mixture of D° and D + , both of which have 
an equal decay rate into 0's. The predicted rates are 

B(B s ^(f)X) = (16.1 ±2.4)% (21) 
B {x>om \B s -> <j>X) = (15.7 ± 2.3)%, (22) 

where we have added in the (1.2±0.4)% from other processes; the rate for 0.05 > x > is taken 
from our Monte Carlo simulation, and amounts to 2.6% of the total yield. 

Solving Eq. (20) and using our procedure for finding the errors by generating Gaussian 
distributions for the independent quantities leads to the probability distributions of fs shown 
in Fig. 9. We measure the B^B S ^* ratio to the total bb quark pair production above the 

1630906-060 
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S 

FIG. 9: Probability distribution of B(T(55) -> B^B^), measured using the inclusive yields of <j) 
mesons, smeared by: (a) statistical error (fitted dots), (b) systematic error (filled triangles) and, (c) 
total error (filled squares). 

four-flavor (u, d, s and c quarks) continuum background at the T(5S) energy as 

B{T(5S) -> B^B { * ] ) = (24.6 ± 2.9™)% . (23) 
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D. Estimate of B(T(5S) -> <f>X) 



We estimate the branching fraction in the 0< x <0.05 interval for T(5S) decays by using 
the ratio of the meson yield in the first x interval to the total. This estimate is obtained 
from a combination of B S B S X and BBX Monte Carlo simulated decays at the T(5S) energy. 
To combine these two types of events, we use the average of our two measurements of the fs 
fraction, 19.3%. Here the fraction in the first bin from B decays is 2.4% very similar to the 
2.6% from B s decays, that makes the error from f s negligible. Adding this estimate to the 
B^ x>om \T{5S) -> <f>X) = (12.9 ± 0.7tU)% as measured in Eq. (14), we find 

B(T(5S) -> cf>X) = (13.8 ± 0.7±?;1)%. (24) 



III. UPDATED MEASUREMENT OF B(T(55) -> B^* ) B S ( * ) ) USING D s MESON 
YIELDS 

Here we update our measurement of f s using D s yields given in Ref. [5] due to changes in 
the scale factors (Eqs. 6 and 7) and therefore to a better estimate of the number of hadronic 
events. We find the ratio 

B(T(5S) -> D S X) , n , 

ff(T(4S) —> D S X) ~ ■ (25) 

We rewrite Eq. (19) for D s mesons rather than for mesons as 



B(T(5S) - D S X) ■ B(D S - 0tt)/2 = f s ■ B(B S - D S X) ■ B(D S - 0tt) (26) 

+ (1 ~ /s) • B(T(AS) ^ ,D S X) • B(D S ^ 0vr) . (27) 

Then using the product of production rates in Ref. [5] updated with the new scale factors, 
our model dependent estimate of B(B S — > D S X) = (92 ±11)% [5] and the newest estimate 

of B(D S -> 0tt) = (4.4 ± 0.6)% [12], we obtain the B^B~ S { * ] ratio to the total bb quark pair 
production above the four-flavor (u, d, s and c quarks) continuum background at the T(5S) 
energy of 

fs = B(T(5S) -> B^B { * ] ) = (16.8 ± 2.6l|I)% (28) 

The probability distribution of this measurement within the statistical, the systematic and 
the total error is shown in Fig. 10. 

This measurement is obtained using a value of (4.4 ± 0.6)% [12] for the D s — > 07T branching 
fraction. Since this branching fraction is not well known, we show in Fig. 11 the behavior of f s 
for different values of B(D S — > tpn) ranging from 3.0 to 5.0%. 

IV. THE T(55) PRODUCTION CROSS SECTION AND B MESON YIELDS 

We measure the cross-section 

a(e + e- -> T(5S)) = (0.301 ± 0.002 ± 0.039) nb , (29) 
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FIG. 10: Probability distributions for B(T(55) -> B^B^), obtained by measuring D s yields, and 
smeared by: (a) statistical error (filled circles), (b) systematic error (filled triangles) and, (c) total 
error (filled squares). 



where the first error is statistical and the second is systematic and due to the error in the 
relative luminosity measurement between T(5S) and continuum. 

Measurement of this cross-section allows us to present the previous CLEO exclusive decay 
results for B mesons at the T(5S) [8] in terms of branching fractions. These are given in 
Table VI. These results are consistent with theoretical expectations [3, 4]. 

The inclusive B meson yield is represented as BBX in the Table. This measurement of 
(58.9 ± 10.0 ± 9.2)% is equal to 1-fs under the assumption that there are no non-66 decays of 
the T(5S). 



V. CONCLUSIONS 

We measure the momentum spectra of mesons from T(4S) and T(5S) decays, and use these 

data to estimate the ratio fs of B^B S to the total bb quark pair production at the T(5S) 
energy as (24.6 ± 2.9l 5 jj)%. We also update our previous measurement of fs using D s yields 
[5], and we find (16.8±2.6jl3;|)%. The central value and the error slightly change due to a better 
determination of the relative luminosity. Furthermore, using our cross-section measurement of 
&5S — (0.301 ± 0.002 ± 0.039) nb for hadron production above 4-flavor continuum at the T(5S) 
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FIG. 11: / s versus B(D S —* c/>7r) (the central line), the statistical error on f s is represented by the inner 
adjacent lines and the total error (the statistical and systematic added in quadrature) is represented 
by the outer two lines. The error on B{D S — > cpn) is not included. 

energy, we measure total B meson production and extract fs = (41.1 ± 10.0 ± 9.2)%. Taking 
a weighted average of all three methods gives fs = (2ljJ)%, where common systematic errors 
have been accounted for. 

Our results as summarized in Table VII, and at the current level of precision, are consistent 
with the previously published phenomenological predictions [3, 4]. 
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TABLE VI: Branching fractions of the possible T(5S) final states either measured directly in this 
paper or estimated using our measurement of the T(5S) total cross section and the cross section and 
upper limits reported in Ref. [8]. 



Different components 
of the T(5S) 



WW 



BBX 
B*B* 

BB* + B*B 
BB 

BB^tt + B^Btt 
BBttit 



Measured 
cross-section (nb) 



Corresponding 
Branching Fraction (%) 



0.177 ±0.030 ±0.016 
0.131 ±0.025 ±0.014 
0.043 ± 0.016 ± 0.006 

< 0.038 (@ 90% CL) 

< 0.055 (@ 90% CL) 

< 0.024 (@ 90% CL) 



(16.8 ± 2.6± 6 3 l)% 
(using D s yields) 

(24.6 ± 2.9^;°)% 
(using <f> yields) 

58.9 ±10.0 ±9.2 
43.6 ±8.3 ±7.2 
14.3 ±5.3 ±2.7 

< 13.8 

< 19.7 

< 8.9 



TABLE VII: Summary of our results. 

Quantity Measurement 
Results from inclusive 4> measurements 
B(T(4S) <f)X) 
B(B -»• <f>X) 
B(T(5S) <j>X) 

B(T(5S) -► 4>X)/B{T{AS) -► 4>X) 

g(T(55) -> WB^) 

Results from inclusive D s measurements 

B(T(5S) -► D S X)/B(T(AS) -► D S X) 2.5 ± 0.2^ 

g(T(55) -> B ( * ] W) (16.8±2.6t|J)% 

Results from cross-section and B measurements 

cr(e + e- -► T(5S)) (0.301 ± 0.002 ± 0.039) nb 

g(T(55) -> B^W) (41.1 ±10.0 ±9.2)% 

/g from combining inclusive 4>, D s and 5 measurements 
g(T(55) B^W) (21±|)% 



(7.06 ±0.09 ±0.60)% 
(3.53 ± 0.05 ± 0.30)% 

(13.8 ± 0.7±?;|)% 
1.9 iO.llg;! 

(24.6 ± 2.9l 1 5 1 ;°)% 
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